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ABSTRACT 


Ice nuclei (IN) significantly affect clouds via supercooled droplets, that in turn 
modulate atmospheric radiation and thus climate change. Since the IN effect is relatively 
strong in stratiform clouds but weak in convective ones, the overall effect depends on the 
ratio of stratiform to convective cloud amount. In this paper, 10 years of TRMM 
(Tropical Rainfall Measuring Mission) satellite data are analyzed to confirm that 
stratiform precipitation fraction increases with increasing latitude, which implies that the 
IN effect is stronger at higher latitudes. 

To quantitatively evaluate the IN effect versus latitude, large-scale forcing data from 
ten field campaigns are used to drive a CRM (cloud-resolving model) to generate long- 
term cloud simulations. As revealed in the simulations, the increase in the net downward 
radiative flux at the TOA (top of the atmosphere) from doubling the current IN 
concentrations is larger at higher latitude, which is attributed to the meridional tendency 
in the stratiform precipitation fraction. 

Surface wanning from doubling the IN concentrations, based on the radiative balance 
of the globe, is compared with that from anthropogenic CCb. It is found that the former 
effect is stronger than the latter in middle and high latitudes but not in the Tropics. 

With regard to the impact of IN on global wanning, there are two factors to consider: 
the radiative effect from increasing the IN concentration and the increase in IN 
concentration itself. The former relies on cloud ensembles and thus varies mainly with 
latitude. In contrast, the latter relies on IN sources (e.g., the land surface distribution) and 
thus varies not only with latitude but also longitude. Global desertification and 
industrialization provide clues on the geographic variation of the increase in IN 
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concentration since pre-industrial times. Thus, their effect on global warming can be 
inferred and then be compared with observations. A general match in geographic and 
seasonal variations between the inferred and observed wanning suggests that IN may 
have contributed positively to global wanning over the past decades, especially in middle 
and high latitudes. 
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1. Introduction 


The global trend in air surface temperature, such as that over the period 1979 to 2005, 
is highly variable from one region to another (IPCC 2007, p. 253). Since carbon dioxide 
(CO 2 ) and other greenhouse gases as well as their effect on global warming are 
horizontally uniform in distribution, it is still unclear how they could cause the observed 
geographical variation in wanning via a feedback mechanism (IPCC 2007; see Section 5 
for more discussion). Thus, other contributions to global warming may exist beyond 
greenhouse gases (e.g., Lindzen 1990; Ramanathan et al. 2007; Prenni et al. 2007). 

Aerosol particles, owing to their direct (Ramanathan et al. 2007; Wild et al. 2007) and 
indirect (NRC, National Research Council 2005; Lohmann and Feichter 2005; Zhang et 
al. 2007) effects on the radiation budget of the Earth, may be responsible for the 
geographical variation in global wanning (Baker 1997; Prenni et al. 2007). Recent cloud- 
resolving model (CRM) simulations have revealed that ice nuclei (IN), a class of aerosol 
particles, can significantly affect cloud ensembles (Phillips et al. 2005; Ekman et al. 
2007; Phillips et al. 2007; Zeng et al. 2008b) and in turn radiation (Zeng et al. 2009, see 
Section 2 for more discussion). An interesting question is whether IN have contributed to 
the geographical variation in global wanning. 

The present paper proposes a connection between increasing IN concentration and 
global wanning. It consists of five sections. In Section 2, TRMM satellite data are 
analyzed to show a meridional variation in stratifonn precipitation fraction and then a 
meridional tendency in the IN effect is proposed. In Section 3, CRM simulations at 
various latitudes are used to confirm and quantify the tendency. In Section 4, the 
geographic variation of IN-induced warming is inferred and the inferred variation 
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compared with observations. Section 5 concludes. 


2. A meridional tendency in the IN effect 

This section as well as the next one addresses meridional and seasonal variations in 
radiation when the current IN concentration is doubled. First, it traces the origin of the IN 
effect downscale to supercooled droplets and then explains why the IN effect relies on the 
ratio between stratiform and convective precipitation rates. Second, it analyzes 10 years 
of TRMM data to quantify the ratio versus latitude, ending with the prediction that there 
exists a meridional tendency in the effects of IN. 

a. Origin of the IN effect 

The IN effect originates with the coexistence of ice crystals and supercooled droplets 
in mixed-phase clouds. Since the saturation water vapor pressure over water is higher 
than that over ice, ice crystals gain mass by vapor deposition at the expense of the 
supercooled droplets that lose mass by evaporation (Wegener 1911). Owing to the 
relatively low ice crystal number concentration but high droplet concentration in clouds, 
many droplets surround each ice crystal. As a result, the droplets evaporate with the 
resulting water vapor being deposited onto crystals to form precipitating particles, which 
is called the Bergeron (1935) precipitation process. 

Consider two cold clouds that have the same structure except that one has IN but the 
other not. Since supercooled droplets exist at temperatures from 0 to -40°C, they are 
thermodynamically unstable. They remain in a liquid state in the cloud without IN. In the 
other cloud with IN, however, some supercooled droplets freeze due to heterogeneous 
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nucleation and others evaporate with the resulting water vapor being deposited onto 


activated IN and ice crystals. This sensitivity of clouds to IN, which is referred to here as 
the Findeisen (1938) colloidal instability, makes it possible for IN to dramatically affect 
clouds and radiation (Jiang et al. 2000; Zeng et al. 2008b, 2009). Since the Bergeron 
process (or Findeisen instability) is one of the dominant precipitation processes in 
stratiform mixed-phase clouds, it leads to IN having a strong effect on stratiform clouds 
(e.g., Jiang et al. 2000). 

b. Effect of IN on stratiform clouds 

Recent long-term CRM simulations revealed that an increase in IN concentration 1 
brings about a decrease in bulk precipitation efficiency (PE) and thus an increase in upper 
tropospheric (UT) cloud ice content (Zeng et al. 2009). This effect of IN on PE can be 
illustrated with a two-layer model. First, consider a single ice crystal that grows by vapor 
deposition and riming in a mixed-phase region. Its depositional growth rate, measured 
with an increase in scale (e.g., equivalent radius), is inversely proportional to its scale, 
whereas the riming growth rate is directly proportional to its scale. Hence, the Bergeron 
process narrows the ice crystal spectrum while riming broadens it. Since the Marshall- 
Palmer distribution of precipitating particles (Pruppacher and Klett 1997) implies the 


1 Owing to the importance of supercooled droplets in the effect of increased IN on 
radiative forcing (see Fig. 11 for a summary), only the increase in IN concentration in the 
mixed-phase region and its effects are modeled and discussed in the present paper and in 
previous ones (Zeng et al. 2008b, 2009). 
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existence of large particles, its widespread occurrence indicates that riming is important 
in mixed-phase precipitation in the atmosphere. 

Second, consider a mixed-phase cloud with a cloud top temperature wanner than - 
40°C. Suppose that all of the ice crystals in the cloud have the same size and the cloud is 
divided vertically into two layers. Given an increase in IN concentration in the upper 
layer, new ice crystals will first form on the activated IN and then begin to approach the 
original ones in size via the Bergeron process, although they remain smaller than the 
original ones. As a result, the precipitating ice content and thus PE increases 
correspondingly (neglecting riming, Graupel riming is quite effective especially when 
graupel particles are large). 

After riming starts to become important, the new crystals cannot approach the original 
ones in size. Instead, their difference in size becomes larger and larger because the riming 
growth rate is directly proportional to crystal scale. Since the new (or smaller) crystals 
grow at the expense of supercooled droplets, cloud water content decreases and therefore 
the growth rate of the original (or larger) crystals slows. In other words, the original 
crystals would decrease in size with increasing IN concentration in the upper layer. Once 
the original crystals fall into the lower layer, the conversion rate of cloud water to 
precipitating ice due to riming would be reduced there as a result of the increased IN in 
the upper layer. 

Suppose that the IN concentration follows the formula of Fletcher (1962). Thus, the 
increase in IN concentration in the lower layer is much smaller than that in the upper one. 
As a result, the change in precipitating ice formation is quite small in the lower layer if 
no precipitating ice particles fall in from the upper layer, but large if precipitating ice 
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crystals fall in. Since the crystal growth rate due to the riming (or Bergeron) process is 
directly (or inversely) proportional to crystal scale, the impact of IN on precipitation 
formation via the riming process is much greater than it is for the Bergeron in the lower 
layer. As a net result, the bulk PE in the mixed-phase cloud decreases with increasing IN 
concentration. 

In summary, with increasing IN concentration, small precipitating crystals (or snow 
content in the CRM) increase in population while large precipitating particles (or graupel 
content in the CRM) decrease. Correspondingly, the Bergeron process becomes stronger 
but the riming becomes weaker, resulting in a net decrease in bulk PE. All of the effects 
of increased IN on snow, graupel and bulk PE are evident in the long-term CRM 
simulations of Zeng et al. (2008b, 2009). 

IN can affect UT ice and radiation via the supercooled droplets in stratiform clouds. 
Consider, for example, a stratiform mixed-phase cloud in the middle troposphere. It is 
usually embedded in a synoptic system (e.g., a cold front; Matejka et al. 1980) with large- 
scale upward motion. Once it dissipates, its residue (e.g., cloud drops) is transported by 
large-scale upward motion into the upper troposphere where it forms ice clouds. Since IN 
can affect stratiform mixed-phase clouds in the middle troposphere, they can modulate 
UT ice and thus radiation on a long timescale. 

c. Effect of IN on convective clouds 

IN can also affect UT ice and radiation via the supercooled droplets in convective 
clouds. They work similarly but in a complicated way, because there are strong updrafts 
and thus high liquid water content in the clouds. Figure 1 displays a schematic cloud 
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cluster (Zipser 1969; Houze 1977) to illustrate how IN affect radiation. Ice crystals and 
supercooled droplets coexist in a mixed-phase region in the middle troposphere. Since IN 
concentration influences ice crystal concentration, it can affect the Bergeron and riming 
processes in the region, which in turn can affect the UT cirrus anvil through rising air 
currents. 

The precipitation formation in convective clouds, just as that in stratiform ones, 
decreases with increasing IN concentration. Since IN affect the Wegener and Bergeron 
processes directly but graupel riming indirectly via the amount of available supercooled 
droplets, the decrease in precipitation formation from increased IN becomes smaller 
when the supercooled droplet content is higher. Consider, for example, an air parcel in a 
convective cloud moving upward starting at the bottom of the lower mixed-phase region. 
Numerous graupel particles that formed in the upper region, as discussed in the preceding 
subsection, fall into the parcel and rime the supercooled droplets. If the IN concentration 
increases, these intruding graupel particles decrease in size and the conversion rate of 
cloud water to precipitating ice (or PE) decreases significantly. 

When the parcel rises slantwise, its temperature falls. As a result, many of the 
increased IN become activated and grow to large ice crystals via the Bergeron and riming 
processes, although they are still smaller than the previous intruding particles. Since the 
parcel has high cloud water content, the new, yet large crystals rime cloud water 
efficiently, partially offsetting the preceding effect of the increased IN on cloud water. As 
a result, the effect of IN on bulk PE in convective clouds is weaker than that in stratiform 
clouds. This difference in the IN effect between convective and stratiform clouds is 
supported by the CRM simulations of Zeng et al. (2008b, 2009). The modeled bulk PE in 
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middle latitudes, for example, decreases more with increasing IN concentration in spring 
than in summer (Zeng et al. 2009) because the stratiform precipitation percentage in 
spring is larger than that in summer. 

d. The ensemble effect of IN on radiation 

The ensemble effect of IN on radiation relies on cloud ensembles as well as individual 
cloud structures. Clouds are usually classified into convective and stratiform (Churchill 
and Houze 1984; Nesbitt et al. 2006). Since the effect of IN on PE is relatively strong in 
stratiform clouds but weak in convective ones, the ensemble effect of IN on radiation 
depends on the ratio between the stratiform and convective cloud amounts, which is 
formulated next. 

Consider a cloud ensemble that consists of two kinds of clouds: stratiform and 
convective. Let N s and N c denote the stratiform and convective cloud amounts, 
respectively. Thus, the increase in downward radiative flux from doubling the IN 
concentration can then be written 

AF, = Af sl N s + Af cl N c , 
or 

A F, = (N s + N c )[Af c , + (A f„ - A f a )F s ], (2.1) 

where the fraction of stratiform clouds F S =N S /(N S +N C ) ; A (/ and A f ci represent the increases 
in the radiative flux from doubling the IN concentration when N s =\ and N c = =1, 
respectively. As shown in (2.1), the fraction of stratiform clouds is important in the 
ensemble effect of IN on radiation because A f si is much larger than Af d in magnitude. 
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Strictly speaking, N s and N c represent the amounts of stratiform and convective mixed- 
phase clouds, respectively. Since mixed-phase clouds usually precipitate efficiently, the 
fraction of stratiform clouds in (2.1) has the same meridional tendency as the fraction of 
stratiform precipitation, although the two fractions differ in quantity to some extent. 
Next, TRMM observations are used to analyze the meridional tendency in stratiform 
precipitation fraction. 

e. Climatology of stratiform precipitation fraction 

Belts of climatic zones imply that cloud properties (e.g., the stratiform precipitation 
fraction) vary mainly with latitude instead of longitude. In this subsection, 10 years of 
TRMM data are analyzed to identify meridional and seasonal variations in stratiform 
precipitation fraction. 

A precipitation radar (PR) was included on the TRMM satellite to provide three- 
dimensional maps of convective and stratiform clouds. In addition, a TRMM microwave 
imager (TMI), a multi-channel passive microwave sensor, was also included to provide 
rainfall information (Simpson et al. 1988). These two instruments have performed 
successfully for over ten years. Their coverage extends from approximately 37°S to 
37°N. 

The latitudinal stratiform precipitation fractions were computed using a gridded (0.5 x 
0.5 degree horizontal resolution) rain product based on TRMM satellite data (i.e., 3G68). 
The product includes rainfall information derived from the TMI and the PR. The TMI 
rain retrieval algorithm (Kummmerow et al. 2001; Olson et al. 2006) uses a Bayesian 
approach to match observed brightness temperatures with a database of CRM-simulated 
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hydrometeor profiles via a radiative transfer model. A combination of techniques using 
polarization information and horizontal textures (i.e., gradients) are used to separate the 
rainfall into its convective and stratiform components (Olson et al. 2001). The PR rain 
retrieval algorithm (Iguchi and Meneghini 1994; Iguchi et al. 2000) accounts for the 
effects of attenuation. Horizontal reflectivity gradients (Steiner et al. 1995) as well as 
bright band information are used to separate the PR rainfall into its convective and 
stratiform components. The stratiform precipitation fraction as a function of latitude is 
calculated by summing the total stratiform rain amount in each 0.5 degree latitude belt 
and dividing by the total rain amount over the same area. Average values were computed 
over a 10-year period ranging from 1 January 1998 through 31 December 2007. 

Figure 2 displays the mean PR- and TMI-derived stratiform precipitation fractions 
against latitude. Since atmospheric stability increases with latitude, the fractions become 
larger at higher latitude and approach 100% over arctic regions (Curry et al. 1996). 
Hence, dashed lines are introduced in the figure to show the meridional tendency beyond 
37°. 

Generally speaking, the PR- and TMI-derived fractions vary similarly with latitude. 
Since the PR beam can penetrate into clouds and thus provide bright band information to 
separate convective and stratiform regions, the PR-derived stratiform precipitation 
fraction makes more sense in Eq. (2.1) than the TMI-derived. The PR-derived fraction, as 
shown in Fig. 2, changes little between 20°S and 20°N but increases significantly with 
latitude in the subtropics and middle latitudes. 

Figure 3 displays the PR- and TMI-derived stratiform precipitation fractions versus 
latitude for all seasons. The fractions are relatively large in winter and spring and 
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relatively small in summer and fall. In particular, the PR-derived fraction in spring is 
close to that in winter in the northern hemisphere and is much larger than that in summer 
and fall in middle latitudes. 

Since the IN effect is relatively strong in stratiform clouds but weak in convective 
ones, Af S i is much larger than A f C j in magnitude. For a given value of ( Af S j-Af C j ), the 
ensemble effect of IN on radiation relies on the stratiform cloud fraction and therefore 
varies with latitude and season via the fraction. Based on Figs. 2 and 3 as well as Eq. 
(2.1), two results are inferred: (1) the ensemble IN effect varies little with latitude in the 
Tropics but significantly in middle and high latitudes and (2) the effect is stronger in 
winter and spring than in summer and fall especially in the northern hemisphere. 

3. Numerical experiments 

In this section, CRM simulations are carried out to confirm and quantify the proposed 
meridional tendency in the IN effect. It consists of four subsections. The first one 
introduces field campaigns as well as their large-scale forcing data. The second describes 
the model structure. The third presents CRM simulations at various latitudes, focusing on 
the modeled meridional variation in stratiform precipitation fraction. The last one 
confirms and quantifies the proposed meridional tendency in the IN effect. 

a. Field campaigns 

Ten field campaigns conducted at various latitudes were used to obtain large-scale 
forcing to drive long-term CRM simulations as well as cloud observations to evaluate the 
modeled cloud ensembles. These field campaigns were scattered over the Tropics and 
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middle latitudes; their geographic locations are shown in Fig. 4. The campaigns were 
chosen primarily to show the insensitivity of radiative forcing (or the net radiative flux at 
the TOA) to IN in the Tropics and the rest to show the sensitivity of radiative forcing to 
IN outside the Tropics. 

Table 1 summarizes the field campaigns and their CRM simulations. In addition to 
the campaigns studied in Zeng et al. (2009), the NAME (North American Monsoon 
Experiment) is added in this study to explore the effects of IN in the subtropics, because 
the stratiform precipitation fraction increases significantly with latitude in the region (see 
Fig. 2). 

NAME was centered at 26°N and 107°W (Johnson et al. 2007). It provided average 
large-scale forcing over the land portion of the NAME Enhanced Budget Array (EBA) 
from 0000 UTC 7 July through 1800 UTC 15 August 2004 (Ciesielski and Johnson 2008; 
Ciesielski et al. 2009). The forcing data were derived from six-hourly observations at five 
sounding sites located along the perimeter of the EBA polygon and included a humidity 
correction (Ciesielski et al. 2009). This forcing data is used to drive 40-day CRM 
simulations to determine the effect of IN in the subtropics. 

In addition to NAME, large-scale forcing data from six tropical field campaigns were 
used to drive CRM simulations to establish the effect of IN in the Tropics. The six 
campaigns include the Tropical Ocean Global Atmosphere (TOGA) Coupled Ocean- 
Atmosphere Response Experiment (COARE) conducted in 1992 and 1993, the 1999 
Kwajalein Experiment (KWAJEX), the 1974 Global Atmospheric Research Program’s 
Atlantic Tropical Experiment (GATE), and the 2006 Tropical Warm Pool - International 
Cloud Experiment (TWP-ICE). The 1998 South China Sea Monsoon Experiment 
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(SCSMEX) that was conducted over the northern (NESA) and southern (SESA) 
enhanced sounding array polygons was also included. 

Three field campaigns were used to provide large-scale forcing to drive CRM 
simulations in middle latiUides. The campaigns were conducted at the Southern Great 
Plains (SGP) site set up by the Atmospheric Radiation Measurement (ARM) Program. 
The three campaigns are referred to here as ARM-SGP-97, 00 and 02, which occurred in 
the summer of 1997, spring of 2000 and summer of 2002, respectively. All of the field 
campaigns except for NAME are summarized in Zeng et al. (2009). 

b. Model description 

The Goddard Cumulus Ensemble model is a CRM (Tao and Simpson 1993; Tao et al. 
2003). Here it is used to simulate the effects of IN on clouds and radiation. The model is 
non-hydrostatic and anelastic. It parameterizes subgrid-scale (turbulent) processes with a 
scheme based on Klemp and Wilhelmson (1978) and Soong and Ogura (1980), and 
incorporates the effects of both dry and moist processes on the generation of subgrid- 
scale kinetic energy. It uses a three-class ice formulation for cloud microphysics based on 
the scheme of Rutledge and Hobbs (1984) with some modification (Lang et al. 2007; 
Zeng et al. 2008b). The model includes the sedimentation of cloud ice (Starr and Cox 
1985) to better simulate clouds in the upper troposphere. It calculates all scalar variables 
(temperature, water vapor, and all hydrometeors) with a positive definite advection 
scheme (Smolarkiewicz and Grabowski 1990). 

The model takes account of both emission and absorption of infrared radiation and 
both scattering and absorption of solar radiation. It divides the entire spectrum of infrared 
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radiation into nine bands to compute infrared fluxes, and thus models the effects of water 
vapor and hydrometeor species on longwave fluxes (Chou et al. 1995, 1999). It also 
divides the entire spectrum of solar radiation into 11 bands (i.e., seven bands in the 
ultraviolet region, one band in the photosynthetically active radiation region, and three 
bands in the near infrared region) to compute solar fluxes, and thus models the effects of 
water species on shortwave fluxes (Chou and Suarez 1999). This radiation scheme is 
summarized in Tao et al. (1996, 2003), and has been used to assess cloud-radiation 
interaction (Tao et al. 1996). 

The model has live prognostic hydrometeor variables: the mixing ratios of cloud 
water, rainwater, cloud ice, snow and graupel. The IN concentration is introduced into the 
parameterization of the Bergeron process as an input factor (Zeng et al. 2008b, 2009). 
Since the IN concentration changes the cloud-ice crystal concentration via heterogeneous 
nucleation, it quickly changes the cloud-ice crystal spectrum because smaller crystals 
grow faster via vapor deposition. Hence, the Bergeron process or the conversion rate of 
cloud ice to snow due to vapor deposition is connected to IN concentration by (Zeng et 
al. 2008b, 2009) 

ma x\2a l (3q i - m I5O p- 1 ^N i )m a I 2 5 -\0'] i (3.1) 

and the Wegener process or the conversion rate of cloud water to cloud ice is connected 
to IN concentration by 

7 77 — [ 3a 2<li + (1 " (3.2) 

{a 2 + 1 )(a 2 + 2) L 

where Nj is the number concentration of active ice nuclei, the mixing ratio of cloud ice, 
a i and a 2 the temperature -dependent parameters in the Bergeron process (Koenig 1971), 
p the air density, the mass of an ice crystal 50 pm in diameter /«/5o = 4.8x10~ g, and the 
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ice particle enhancement factor due to a riming-splintering mechanism (Hallett and 
Mossop 1974) and others //= 1 .2 in mid-latitudinal clouds. Expressions (3.1) and (3.2) 
follow the cgs (centimeter-gram-second) system, and the two conversion rates have units 
ofgg' 1 s' 1 . 

Since the ice crystal concentration is much lower than water droplet concentration in 
the mixed-phase region, the air is almost saturated with respect to water (e.g., Korolev 
and Mazin 2003). Hence, the diffiisional growth of ice crystals is independent of water 
droplet concentration. Therefore, (3.1) and (3.2) are independent of cloud water content. 

The model employs the formula of Fletcher (1962) to compute the active IN 
concentration in clouds as a function of air temperature T or 

N i =n 0 Qxp[fi(T 0 -T)\ (3.3) 

where n 0 and p are constant, and 7o=273.16 K. Substituting (3.3) into (3.1) and (3.2) 
yields the two expressions for the Bergeron and Wegener processes, respectively. 
Obviously, the expressions vary with firio and /?, where /i and uq are not separated. 

To explore the effect of IN on clouds and radiation, the model is run with differing 
amounts (i.e., categories) of IN concentration (see Table 2) for a given large-scale 
forcing. The ice crystal number concentration, generally speaking, increases as the IN 
category changes from I to VI. 

Since the model represents clouds explicitly (Moncrieff et al. 1997; Randall et al. 
2003), it is used to simulate clouds with prescribed large-scale forcing derived from 
observational data. The simulations follow the same setup as has been used in previous 
studies (e.g., Grabowski et al. 1998; Petch and Gray 2001; Blossey et al. 2007; Zeng et al. 
2007, 2008b, 2009) of this type. All of the simulations carried out are three-dimensional 
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(3D), using a 1 km horizontal resolution and a vertical resolution that ranges from 42.5 m 
at the bottom to 1 km at the top. The model uses 256x256x41 gridpoints and a time step 
of 6 seconds for integration. 

c. CRM simulations at various latitudes 

CRM simulations for the ten field campaigns are carried out to address the IN effect 
versus latitude. Different IN categories (see Table 1) are tested for each campaign. Most 
of the simulations are detailed in Zeng et al. (2009). Here, four additional CRM 
simulations for NAME are introduced to bridge the gap between the Tropics and middle 
latitudes. The four simulations follow the same model setup as those in Zeng et al. (2009) 
except for the following details. The simulations are driven with the large-scale forcing 
derived from the NAME observations (Ciesielski et al. 2009), starting at 0000 UTC 7 
July 2004 and lasting for 40 days. Since it is difficult to get accurate mean surface fluxes 
over a land region for long-term CRM simulations (Zeng et al. 2007), the simulations 
relax the temperature and humidity at the lowest model level to the observed values with 
a timescale of ten minutes. This relaxing timescale is so short that the model results are 
not sensitive to the imposed surface fluxes. 

The four simulations use IN categories II to V (see Table 2), respectively. The 
modeled precipitation rate and precipitable water are sensitive to IN concentration (figure 
omitted) just as in the other CRM simulations (Zeng et al. 2009). To compare 
precipitation formation in the stratiform mixed-phase region with that in the convective 
one, a contoured frequency with altitude diagram (CFAD; Yuter and Houze 1995) of 
vertical velocity in the mixed-phase region is computed for the NAME simulation that 
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uses IN category IV (figure omitted), showing that small vertical velocities are far more 
prevalent than large ones in the mixed-phase region. 

Since the variables q / and q c q g are directly proportional to the precipitation formation 
rates due to the Bergeron and riming processes, respectively. Fig. 5 displays their 
averages as functions of pressure and vertical velocity, where q c and q g are the mixing 
ratios of cloud ice, cloud water and graupel, respectively. The values in the figure are 
computed by averaging q and qcq g horizontally in the same mixed-phase region. As 
shown in Fig. 5, graupel riming is relatively important in the regions with large vertical 
velocity (or convective precipitation). In contrast, the Bergeron process is important in 
the regions with small vertical velocity (or stratiform precipitation) because the frequency 
of low vertical velocity is relatively high. 

The NAME simulations are consistent with the others in duplicating the shift of 
precipitation from convective to stratiform with increasing latitude. Figure 6 displays the 
averages of q , and qcq g as functions of pressure and vertical velocity for two CRM 
simulations: one over a tropical region (i.e., the GATE simulation with IN category V 2 * ) 
and the other over a mid-latitudinal region (i.e., the ARM-SGP-00 simulation with IN 
category II). The tropical simulation, in contrast to the mid-latitudinal one, has strong 


2 The proper IN category (or /W 7 ) for each campaign is chosen so that its corresponding 

CRM simulations are close to the field observations (Zeng et al. 2009; also see section 

3.e). In the present paper, N t is not distinguished from //. In fact, //decreases significantly 
with increasing latitude. Its magnitude in tropical clouds is about 10M0 4 times that in 
mid-latitudinal clouds, which is discussed elsewhere. 
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graupel riming especially at a pressure of 420 hPa and vertical velocity of 25 m/s, which 
indicates that graupel riming is more important in tropical precipitation. 

Contrasting Figs. 5 and 6 shows that graupel riming is most important in the GATE 
simulation, less so for NAME, and the least for ARM-SGP-00, which makes physical 
sense because, latitudinally, NAME is located between GATE and ARM-SGP. Table 3 
summarizes the modeled stratiform precipitation percentages in GATE, NAME and 
ARM-SGP-00. It shows that the modeled stratiform precipitation fraction becomes larger 
at higher latitude, which agrees with the TRMM observations (or Fig. 2). Hence, it is 
expected that CRM simulations at various latitudes can capture the proposed meridional 
variation in the IN effect. 

d. The IN effect at different latitudes 

IN affect clouds and radiation via supercooled droplets (see Sections 2.b and 2.c). An 
increase in IN concentration brings about a decrease in bulk PE and thus an increase in 
UT cloud ice content (Zeng et al. 2009). Figure 7 summarizes the modeled effects of IN 
on UT cloud ice content and radiation at various latitudes. Since the mixed-phase region 
is observed at temperatures between 0 and -40°C (and most often at the wanner end of 
the range; e.g., Pruppacher and Klett 1997; DeMott et al. 2003), the figure uses the IN 
concentration at -10°C as an index to represent the concentration of active IN in clouds. 

Figure 7 displays the modeled UT cloud ice content against the IN concentration at - 
10°C over the ten field campaigns. The content increases with IN concentration. The 
figure also displays the modeled downward net radiative flux at the TOA versus the IN 
concentration. The flux increases first and then decreases considerably with increasing IN 
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concentration in mid-latitudes. In contrast, the flux changes little with IN concentration in 
the Tropics. 

This differing IN effect can be understood with a simple model. Consider an optically- 
thin UT cirrus cloud in a clear sky. The cloud is so thin that it functions like a subvisual 
cirrus cloud. Thus, its solar reflection is quite small while its infrared absorption is 
relatively large. When its ice content increases, its solar reflection increases slower than 
its infrared absorption. As a result, the downward net radiative flux at the TOA increases 
correspondingly. This sensitivity of the TOA radiative flux to UT ice content explains the 
mid-latitudinal phenomenon that the TOA flux, when the IN concentration or UT ice 
content is low, increases with increasing IN concentration. 

If the ice content continues to increase until the infrared absorption is close to one but 
the solar reflection is still far below one, the ice content would be high. If the content 
increases further, the infrared absorption increases little because of its ceiling, whereas 
the solar reflection can increase a lot. As a result, the downward net radiative flux at the 
TOA decreases correspondingly, which explains the mid-latitudinal phenomenon that the 
TOA flux, when the IN concentration or UT ice content is high, decreases with increasing 
IN concentration. 

In fact, the TOA radiative flux usually changes more moderately compared to those in 
the preceding two extreme cases. When the IN concentration or UT ice content is 
relatively low, the TOA flux increases with increasing IN concentration. When it is high, 
the TOA flux decreases. Of course, when it is moderate, the TOA flux changes little, 
which explains the insensitivity of the TOA radiative flux to IN concentration in the 
Tropics. 
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e. Meridional variation in radiative forcing from doubling the current IN concentration 

IN concentration varies significantly from one region to another (e.g., Fletcher 1962; 
Pruppacher and Klett 1997; DeMott et al. 2003). Given in-situ active IN concentrations 
over the field campaigns, the IN effect from doubling the current IN concentration can be 
determined based on Fig. 7. However, no direct IN observations are available from the 
field campaigns. To diagnose an in-situ IN concentration over a region, the model was 
assigned different IN concentrations for different simulations. If a cloud simulation can 
be made to match the associated observations, then its assigned IN concentration should 
correspond to the in-situ IN concentration over that region. 

To determine the average IN concentration in mid-latitudinal clouds, ARM-SGP 
simulations using different IN concentrations are compared with field observations first. 
Figure 8 displays the biases in the modeled upward infrared flux at the TOA, precipitable 
water and surface precipitation rate against the IN concentration at -10°C. Corresponding 
zero biases indicate that the in-situ active IN concentration is ~3xl0' 2 m’ 3 at -10°C in 
mid-latitudinal clouds. This low IN concentration is consistent with previous 
observations' 1 (e.g., Rutledge and Hobbs 1984) because in many mid- latitudinal 

1 A typical IN concentration measured outside of mid-latitudinal clouds is 1 to 10 L’ 1 at 
temperatures from -17 to -50°C (DeMott et al. 2003; Richardson et al. 2007). It is much 
higher than the diagnosed IN concentration of -3x10° L’ 1 at -10°C in mid-latitudinal 
clouds, even after being adjusted to the same temperature. Since the measured ambient 
IN concentration cannot explain the lack of ice crystals in many mid-latitude continental 
clouds with cloud top temperatures below 0°C (e.g., Paluch 1979; Blyth et al. 1988), it 
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continental cumulus clouds, for example, no ice crystals have been found even with a 
cloud top temperature of -20°C (e.g., Paluch 1979; Blyth et al. 1988). 

Since the diagnosed IN concentration is between category I and II, the ARM-SGP 
simulations with categories I and II are used to determine the IN effect from doubling the 
current IN concentration in middle latitudes. The increase in a variable F (e.g., a radiative 
flux) from doubling the IN concentration is computed as 

p n _ p 1 p u _ p 1 

ln(JVf/jV') ln2 = log, [(>Arf )/(,<)]’ (3 ' 4) 

where a superscript indicates an IN category used in a CRM simulation. Using (3.4), the 
increases in the upward solar and infrared fluxes at the TO A from doubling the IN 
concentration are first obtained from ARM-SGP-97, -00 and -02 and then displayed 
against latitude in Fig. 9. Similarly, the increases in the downward net radiative flux at 
the TOA from doubling the IN concentration are displayed in Fig. 10. 

The preceding procedure for determining the current IN effect is applied to cases over 
the Tropics and subtropics in turn. The modeled precipitation rate and precipitable water 
are compared with the observed over the tropical campaigns, showing that the 
simulations with category V best mimic the observations (Zeng et al. 2009). Recently, 
TRMM PR observations over KWAJEX (Matsui et al. 2009) and radar observations over 
TWP-ICE (to be reported later) were compared with modeled cloud ensembles, 
supporting the conclusion of Zeng et al. (2009) that the tropical simulations with IN 

must be different from the active IN concentration in clouds. Physically this makes sense 
because cloudy and ambient parcels move to their positions along different trajectories 
and therefore should have different rates (or “amounts”) of nucleation. 
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category V resemble observations (see footnote 2 for the role of //). Hence, the CRM 
simulations with IN categories IV and V over the tropical campaigns (i.e., 
SCSMEX/NESA, TWP-ICE, GATE, KWAJEX, SCSMEX/SESA and TOGA-COARE) 
are used to determine the increases in the radiative fluxes from doubling the current IN 
concentrations in the Tropics. Since NAME is located between the ARM-SGP and 
tropical campaigns, the two NAME simulations with IN categories II and IV are used to 
determine the increases in the radiative fluxes from doubling the IN concentration in the 
subtropics. 

Figures 9 and 10 summarize the effect of doubled IN on radiation at different 
latitudes. The upward solar and infrared fluxes at the TOA increase and decrease with 
doubled IN concentrations, respectively, which is attributed to the increase in UT ice 
content with increasing IN concentration. Moreover, the increase in the solar flux is close 
to the decrease in the infrared one in the Tropics, whereas the former is much smaller 
than the latter in middle latitudes. The differing response between the Tropics and middle 
latitudes makes physical sense because the in-situ UT ice content is much higher in the 
Tropics than in middle latitudes (see Fig. 7). 

Figure 10 shows the increase in radiative forcing (or the downward net radiative flux at 
the TOA) from doubling the current IN concentration against latitude, which is obtained 
from the CRM simulations over the ten field campaigns. Based on the meridional 
tendency in stratiform precipitation fraction (see Fig. 2), two lines are introduced to fit 
the simulated increase in spring and summer, respectively. As shown in the figure, the 
current downward net radiative flux at the TOA increases significantly with IN 
concentration in mid-latitudes but not in the Tropics. Furthermore, the net flux is more 
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sensitive to IN concentration in spring than in summer. This differing IN effect in 
different seasons is consistent with that at different latitudes, because stratiform 
precipitation fraction is larger in spring than in summer (see Figs. 2 and 3). 

4. Increasing IN concentration and global warming 

a. Meridional tendency of the IN-induced warming 

IN affect global wanning via the net radiative flux at the TOA. Consider an increase 
in IN concentration. It results in a A F increase in global mean (downward) net radiative 

flux at the TOA via clouds. It also results in an increase in the equivalent blackbody (or 

surface) temperature of the globe from T to T+ AT . Since the Earth maintains radiative 
equilibrium on a climatic scale 4 , 

AF = cr(T + AT) 4 - crT (4.1) 

where a is the Stefan-Boltzmann constant, which yields 

AT=XAF (4.2) 

where X = ( 4<t 7’ 3 ) -1 . The value of X is around 0.25 K (W m' 2 )' 1 (Chamey 1979). 


4 The effective thermal capacity of the Earth’s surface is not important for global 
warming. Consider a linear model C dAT/dt = A F - Act T AT, where t is time and C is the 
effective thermal capacity of the Earth’s surface that wanning influences. For a given 
forcing A F ~ sin {cot) where co is angular frequency, the forced mode is 

AT- |4<rr 3 sm(a>t)-Ca>cos(cot) Jc ©) 2 + (AcrT*) 2 J . When 4erT /Cco» 1, the forced 
mode degenerates into Eq. (4.1), and therefore the thermal capacity becomes negligible. 
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Since radiative equilibrium approximately describes a regional thermal balance 
(Manabe and Strickler 1964; Emanuel et al. 1994), Eq. (4.2) can be used to estimate the 
temperature change for a given regional increase in IN concentration. Supposing that 
T=0.25 K (W m' 2 )’ 1 (Chamey 1979), the vertical axis of Fig. 10 also scales to the increase 
in surface temperature from doubling the IN concentration. 

To compare the IN-induced warming with the CCE-induced one, Figure 10 displays 
the increase in the net radiative flux (or surface temperature) when the atmospheric C0 2 
concentration is increased from a pre-industrial value of about 280 ppm to 379 ppm in 
2005 (IPCC 2007, p. 25). Clearly, the surface wanning from doubling the IN 
concentration is much stronger than that from the anthropogenic CO 2 in the subtropics 
and middle latitudes, and a stronger warming is expected in high latitudes. 

b. Geographic distribution of IN sources 

The IN-induced wanning varies geographically because it depends on two factors: the 
radiative effect from doubling the IN concentration and the actual increase in IN 
concentration. Since the C0 2 -induced warming is horizontally uniform, the IN-induced 
wanning can be identified from global wanning if it varies in space and time similar to 
what is observed. Next, the geographic and seasonal variations in IN-induced warming 
are inferred and contrasted with observations. 

IN (such as dust, biological and soot particles) are usually generated by natural and 
human activity over continents (Pmppacher and Klett 1997, DeMott et al. 2003; Mohler 
et al. 2007; Phillips et al. 2008). Assuming that a significant increase in IN production is 
due to global desertification (Schlesinger 1990), a large increase in IN concentration 
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would arise over continents but not over oceans, which could then lead to a strong 
wanning response over continents but not over oceans. Because the net radiative flux at 
the TO A is insensitive to IN concentration in the Tropics, the warming is weaker over the 
equator but stronger away from the equator. This saddle-shaped warming pattern 
resembles the observed, such as that over the African continent and the American and 
Asian-Australian regions (IPCC 2007, p. 37). As a special case, IN generated in Asia and 
Australia are usually transported by westerly winds over the Pacific such that their 
increase brings about a saddle-shaped warming pattern over the western and central 
Pacific (IPCC 2007, p. 37). 

In addition, localized areas of intense desertification can lead to strong local increases 
in IN production and correspondingly to localized areas of strong wanning. Deforestation 
and desertification near the borders of the Sahara and Gobi deserts, for example, have 
increased IN in those regions, which is well correlated with the strong warming that’s 
been observed over northern Africa and northeastern China (IPCC 2007, p. 37), 
respectively. The central Sahara and western Australia, due to their sparse flora, provide 
little room for further desertification (or no increase in IN production), which results in 
relatively weak wanning in those regions (IPCC 2007, p. 37). 

Besides global desertification, industrialization (Ramanathan et al. 2007) has 
increased IN concentrations (e.g., Plaude et al. 2007) that can contribute to global 
wanning. Consider, for example, a horizontally uniform increase in IN concentration due 
to industrialization. Owing to the meridional tendency in the IN effect, it can bring about 
a stronger warming at higher latitudes, which matches the observation that the strongest 
warming has occurred in the northern arctic region over the past 100 years (IPCC 2007, 
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p. 37). Similarly, heavy industrial IN production over Europe corresponds to the observed 
strong warming there from 1979 to 2005 (IPCC 2007, p. 37). 

The IN-induced warming has a seasonal variation, too. Since the current net radiative 
flux in mid-latitudes would increase more with IN concentration in spring than in 
summer, a similar increase in IN concentration can bring about a stronger wanning in 
winter and spring but a weaker one in summer and fall. This IN effect resembles the 
observations from the last 30 years that the greatest wanning has occuned in winter and 
spring (IPCC 2007, p. 37). All of the preceding scenarios wherein strong warming 
coincides with increased IN production suggest that IN may have contributed to global 
warming especially in middle and high latitudes. 

5. Conclusions and discussion 

The total surface temperature increase from the period 1850 to 1899 to the period 
2001 to 2005 is 0.76°C (IPCC 2007, p. 36), and the contribution of anthropogenic CO 2 to 
the net radiative flux at the TOA is 1 .66 W m'~ in 2005 (IPCC 2007, p. 32). Given T=0.25 
K (W m' 2 )’ 1 in (4.2) (Chamey 1979), the increase in CO 2 over the past 150 years brings 
about a temperature increase of 0.4 1°C or about half of the measured increase. Hence, 
other candidates may exist that could explain the difference. 

Global warming is not horizontally uniform. IPCC (2007, p. 37) summarized its 
geographic and seasonal variations as follows: 

“ Surface temperatures over land regions have warmed at a faster rate than over 
the oceans in both hemispheres. ... The warming in the last SO years is 
widespread over the globe, and is greatest at higher northern latitudes. The 
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greatest warming has occurred in the nor them -hemispheric winter (DJF) and 
spring (MAM). ” 

Based on these variations as well as their details (e.g., those over the period 1979 to 
2005; IPCC 2007, p. 37), the IN-induced warming can be identified in the observed total 
wanning because the effect of C0 2 and other greenhouse gases on global wanning is 
geographically uniform. In reference to the following conclusions, increasing IN 
concentration is proposed as a candidate for global wanning. 

• IN affect clouds and radiation significantly via supercooled droplets (see Fig. 11 for a 
summary). Since the IN effect is relatively strong in stratiform clouds but weak in 
convective ones, the ensemble effect of IN on radiation depends on the ratio between 
stratifonn and convective cloud amount, which is formulated in (2.1). 

• Ten years of TRMM data reveal that the stratifonn precipitation fraction increases 
significantly with latitude beyond the Tropics. The fraction is larger in winter and spring 
than in summer and fall in middle latitudes, especially in the northern hemisphere. 

• CRM simulations at various latitudes show that the increase in the net TO A radiative 
flux from doubling the IN concentration becomes stronger at higher latitude. Moreover, 
the increase is larger in spring than in summer in middle latitudes. These modeling results 
confirm and quantify the meridional tendency of the IN effect derived from the preceding 
two conclusions. 

• It is well-known that IN come from land regions instead of oceans. Thus, global 
desertification and industrialization over the past decades could bring about a large 
increase in IN concentration over land regions but not over oceans, although no 
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quantitative description is available on the geographic variation of increased IN right 


now. 

Superimposing the last two results yields geographic and seasonal variations of the 
IN-induced wanning (see Section 4.b for details). A general match in the variations 
between the infened and observed warming suggests that IN may have contributed 
positively to global wanning especially in middle and high latitudes. With more and more 
information on geographic variations in IN becoming available in the future, global cloud 
modeling systems (e.g., the multi-scale modeling framework; Grabowski 2001) can be 
used to further evaluate the proposed connection between increased IN and global 
wanning. 

On the other hand, regional CRM simulations can still be used to quantify the effect of 
IN on warming over specific geographic regions such as the Arctic. In this study, the 
effects of IN on warming in the Tropics and middle latitudes are extended into polar 
regions based on the meridional variation in stratifonn precipitation fraction. Due to the 
poor performance of current CRMs in modeling arctic clouds and precipitation (e.g., 
clear-sky ice precipitation), no effect of IN on wanning has been quantified in the Arctic. 
Recently, Zeng (2008a) proposed a new precipitation mechanism that works in thin ice 
clouds and even clear skies. With a proper representation of the mechanism, it is possible 
to use a CRM to quantify the effect in the Arctic. Since the new mechanism is not as 
sensitive to IN concentration as the Bergeron and riming processes, it cannot violate the 
proposed connection between increased IN and warming in the Arctic, although the 
discussion is still open. 
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Furthermore, IN observations and modeling still remain as a challenge. In the present 
study, the contribution of increased IN to global warming is estimated based on the 
consistency in geographic variation between the IN-induced and observed warming. 
Since IN concentration varies highly from one region to another, the estimated 
contribution is contaminated with uncertainties. Hence, more and better IN measurements 
are needed as essential information to quantitatively assess the contribution of IN to 
global wanning. 
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Caption 


Figure 1 Schematic of a cloud cluster with an anvil. A cirrus anvil extends laterally 
outward from the precipitating core of the cluster and modulates solar and infrared 
radiation. Ice crystals and water droplets, denoted with stars and circles respectively, 
coexist in the middle troposphere and thus provide an enviromnent for the Bergeron and 
riming processes. Updrafts transport surviving supercooled droplets into the upper 
troposphere where they then freeze. With increasing IN concentration, the surviving 
supercooled droplets in the middle troposphere and thus the ice crystals in the upper 
troposphere increase correspondingly (see Sections 2.b and 2.c for details). 

Figure 2 Meridional variations in stratiform precipitation fraction derived from 10 years 
of TRMM PR and TMI data. Dashed lines are estimated to represent the tendency beyond 
37°. 


Figure 3 Meridional variations in stratiform precipitation fraction derived from 10 years 
of TRMM PR (top) and TMI data (bottom) for four seasons: northern-hemispheric winter 
(thick solid), spring (thick dashed), summer (thin solid) and fall (thin dashed line). 

Figure 4 Geographic locations of the ten field campaigns from which data are used to 
drive and evaluate the present CRM simulations. 

Figure 5 Mean mixing ratio of cloud ice q t (left) and q^g, the product of the mixing 
ratios of cloud water and graupel (right), as functions of pressure and vertical velocity w 
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in the mixed-phase region of the NAME simulation with IN category IV. The units of g, 
and q c q g are g/kg and 1 0' (g/kg)“, respectively. 

Figure 6 As in Fig. 5, but for the GATE simulation with IN category V (top) and the 
ARM-SGP-00 simulation with IN category II (bottom). 

Figure 7 Upper-tropospheric ice content and TOA radiative flux versus IN concentration 
over various geographic regions. Modeled cloud ice content above 7.4 km (top) and net 
downward radiative flux at the TOA (bottom) vary with IN concentration, which are 
obtained from the CRM simulations over the ten field campaigns. One line corresponds 
to one field campaign. Red and blue lines display the mid-latitudinal results in spring and 
summer, respectively; green and black lines display sub-tropical and tropical results, 
respectively. The horizontal axis represents the IN concentration calculated with the 
Fletcher formula at a temperature of -10°C. Red labels correspond to red lines. 

Figure 8 Model biases versus assigned IN concentration in mid-latitudes. Biases in 
modeled upward infrared flux at the TOA (thin lines in the lower panel), precipitable 
water (thick lines) and surface precipitation rate (thin lines in the upper panel) vary with 
the assigned IN concentration at -10°C. A zero bias corresponds to an in-situ active IN 
concentration. Red and blue lines represent the results in spring and summer, 
respectively. 
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Figure 9 Increase in the upward solar flux (black) and decrease in the upward infrared 
flux (red) at the TO A from doubling the IN concentration versus latitude. All of the 
results are obtained from the CRM simulations for the ten field campaigns. 

Figure 10 Increase in the net radiative flux at the TOA from doubling the current IN 
concentration versus latitude. All of the results are obtained from the CRM simulations 
over the ten field campaigns. Thick and thin solid lines are introduced to fit the results to 
spring and summer, respectively, based on Figs. 2 and 3. The dashed line represents the 
increase in the net radiative flux when the atmospheric CO 2 concentration is increased 
from a pre-industrial value of about 280 ppm to 379 ppm in 2005. The vertical axis also 
scales to the increase in surface temperature from doubling the IN concentration when 

M25K(Wm' 2 )''. 

Figure 11 A flowchart showing the effect of global desertification and industrialization 
on warming via IN and supercooled droplets, where IN affect the Wegener and Bergeron 
processes directly and graupel riming indirectly via the amount of available supercooled 
droplets. All the processes in the blue area are explicitly modeled in the present CRM 
simulations. 
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Table 1 Field campaigns and corresponding CRM simulations 


Field Campaign 

Latitude 

Starting Date 

Modeling Days 

IN Category 

ARM-SGP-97 

37°N 

18 Jun 1997 

29 

I, II, IV, V 

ARM-SGP-00 

1 Mar 2000 

20 

I, II, IV, V 

ARM-SGP-02 

25 May 2002 

20 

I, II, III, IV, V 

NAME 

26°N 

7 July 2004 

40 

II, III, IV, V 

SCSMEX/NESA 

21°N 

6 May 1998 

44 

IV, V, VI 

TWP-ICE 

12°S 

19 Jan 2006 

24 

IV, V 

GATE 

9°N 

1 Sep 1974 

18 

IV, V, VI 

KWAJEX 

9°N 

24 Jul 1999 

52 

IV, V 

SCSMEX/SESA 

5°N 

1 May 1998 

60 

IV, V 

TOGA-COARE 

2°S 

1 Nov 1992 

61 

IV, V 
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Table 2 IN concentration categories used in numerical experiments 


IN Category 

jUtio (cm' 3 ) 

P 

I 

6.0xl0’ 10 

0.3 

II 

1.2x1 O' 9 

0.4 

III 

4.8xl0' 9 

0.45 

IV 

1.2x1 O’ 8 

0.5 

V 

1.2x1 O' 6 

0.6 

VI 

1.2x1 O' 5 

0.7 
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Table 3 Modeled meridional variation in stratiform precipitation percentage 


Experiment 

GATE 

NAME 

ARM-SGP-00 

Latitude 

9°N 

26°N 

37°N 

Stratiform Precipitation 

35.0 

44.7 

53.5 

Percentage (%) 
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Figure 1 Schematic of a cloud cluster with an anvil. A cirrus anvil extends laterally 
outward from the precipitating core of the cluster and modulates solar and infrared 
radiation. Ice crystals and water droplets, denoted with stars and circles respectively, 
coexist in the middle troposphere and thus provide an environment for the Bergeron and 
riming processes. Updrafts transport surviving supercooled droplets into the upper 
troposphere where they then freeze. With increasing IN concentration, the surviving 
supercooled droplets in the middle troposphere and thus the ice crystals in the upper 
troposphere increase correspondingly (see Sections 2.b and 2.c for details). 
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Figure 2 Meridional variations in stratiform precipitation fraction derived from 10 years 
of TRMM PR and TMI data. Dashed lines are estimated to represent the tendency beyond 
37°. 
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Figure 3 Meridional variations in stratiform precipitation fraction derived from 10 years 
of TRMM PR (top) and TMI data (bottom) for four seasons: northern-hemispheric winter 
(thick solid), spring (thick dashed), summer (thin solid) and fall (thin dashed line). 
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ARM-SGP 



FIELD CAMPAIGNS 


Figure 4 Geographic locations of the ten field campaigns from which data are used to 
drive and evaluate the present CRM simulations. 
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q NAME (3D) q^, NAME (3D) 



Figure 5 Mean mixing ratio of cloud ice q t (left) and q^g, the product of the mixing 
ratios of cloud water and graupel (right), as functions of pressure and vertical velocity w 
in the mixed-phase region of the NAME simulation with IN category IV. The units of qj 
and q c q g are g/kg and 10' 1 (g/kg) 2 , respectively. 
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Pressure (hpa) 


GATE (3D) 


GATE (3D) 


Qi 


q c q g 



w (m/s) w (m/s) 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


Figure 6 As in Fig. 5, but for the GATE simulation with IN category V (top) and the 


ARM-SGP-00 simulation with IN category II (bottom). 
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Figure 7 Upper-tropospheric ice content and TOA radiative flux versus IN concentration 
over various geographic regions. Modeled cloud ice content above 7.4 km (top) and net 
downward radiative flux at the TOA (bottom) vary with IN concentration, which are 
obtained from the CRM simulations over the ten field campaigns. One line corresponds 
to one field campaign. Red and blue lines display the mid-latitudinal results in spring and 
summer, respectively; green and black lines display sub-tropical and tropical results, 
respectively. The horizontal axis represents the IN concentration calculated with the 
Fletcher formula at a temperature of -10°C. Red labels correspond to red lines. 
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Figure 8 Model biases versus assigned IN concentration in mid-latitudes. Biases in 
modeled upward infrared flux at the TOA (thin lines in the lower panel), precipitable 
water (thick lines) and surface precipitation rate (thin lines in the upper panel) vary with 
the assigned IN concentration at -10°C. A zero bias corresponds to an in -situ active IN 
concentration. Red and blue lines represent the results in spring and summer, 
respectively. 
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Figure 9 Increase in the upward solar flux (black) and decrease in the upward infrared 


flux (red) at the TO A from doubling the IN concentration versus latitude. All of the 


results are obtained from the CRM simulations for the ten field campaigns. 
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Decrease in Upward Infrared Flux at TOA (W/m 2 ) 



Figure 10 Increase in the net radiative flux at the TOA from doubling the current IN 
concentration versus latitude. All of the results are obtained from the CRM simulations 
over the ten field campaigns. Thick and thin solid lines are introduced to fit the results to 
spring and summer, respectively, based on Figs. 2 and 3. The dashed line represents the 
increase in the net radiative flux when the atmospheric CO 2 concentration is increased 
from a pre-industrial value of about 280 ppm to 379 ppm in 2005. The vertical axis also 
scales to the increase in surface temperature from doubling the IN concentration when 
/t=0.25 K (W m" 2 )' 1 . 
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Figure 11 A flowchart showing the effect of global desertification and industrialization 
on wanning via IN and supercooled droplets, where IN affect the Wegener and Bergeron 
processes directly and graupel riming indirectly via the amount of available supercooled 
droplets. All the processes in the blue area are explicitly modeled in the present CRM 
simulations. 


57 









